A Cosmological Test of Standard Gravity by Weak Lensing 
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The large scale structure formation in the standard Einstein gravity at later time is uniquely 
determined by the expansion history (2) measured by the geometrical factor. A possible departure 
from the expected standard structure formation is detectable by growth structure test. We design 
new deviation parameters to test standard gravity at large scales and estimate the expected errors 
by weak lensing experiment. The effective Newtonian constant Gjv at the present time can be 
constrained around 10% accuracy by measuring the difference curvature perturbations ($ — *t)/2 
within 2% accuracy. 
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Introduction The gravity known to us by direct measure- 
ment ranges just from mm to the solar system scales. 
Considering the cosmological gravitational effects beyond 
those scales, we are still free to imagine a possible broken 
standard gravity. The observed cosmic acceleration []| 
can be the evidence of a departure from standard gravity 
at large scales. With the presence of standard gravity, 
cosmic acceleration can be explained by introducing an 
exotic energy component like dark energy (DE). With the 
absence of standard gravity, cosmic acceleration can be 
generated by modified gravity (MG). The extra dimen- 
sion can be introduced to generate cosmic acceleration 
on a 4D brane world (DGP) |^ ,^J. Or we can modify 
the Ricci scalar R with correction terms which are pro- 
portional to the inverse of i? 0. Cosmic acceleration 
is probed by the gravitational effects like the geometri- 
cal factor (GEF) and the growth factor (GRF). By using 
the simultaneous fitting GEF and GRF, we are able to 
distinguish DE and MG 

In this letter, we pay attention to test standard gravity 
at large scales. In standard gravity paradigm, observable 
GRF by weak lensing is uniquely determined by mea- 
sured GEF precisely. The expected fixed evolution of 
GRF leads us to think the way to detect the deviation 
from standard gravity. Inspired by the DGP type MG 
models, we design deviation parameters to probe that de- 
parture. We add the deviation parameters to the usual 
DE cosmological parameters. If we find the fixed GRF 
expected marginally by GEF measurement, then we can 
claim the evidence of Einstein gravity at largest scales. 
But if we find any significant departure from the expected 
DE GRF, then we can argue that standard gravity is bro- 
ken down at large scales (not necessarily excluding DE). 

We estimate the expected deviation parameter errors 
by using cosmic shear experiments. The shapes of galax- 
ies are distorted by the intervening lensing potential 
formed by GRF and projected by GEF. Thus it is a best 
combined test of GEF and GRF. The future weak lens- 
ing experiment like G27r puts a tight constraint on GRF 



within 2% accuracy which leads to the precise test of 
standard gravity at large scales. We find that the effec- 
tive Newtonian constant Gj\j can be constrained around 
10% level of precision. 

GEF test We fit DGP type MG models to find out if 
GEF test alone can test standard gravity broken down. 
The combined test of SN and CMB is considered here. 
We use the SN gold set and the angular diameter 
distance D{ziss) to the last scattering surface measured 
by WMAP .lA]- We fix the distance to recombination 
at ziss = IO88I2 calculated by the acoustic peak scale 
£a = 299±2 and its length calibration through the matter 
density fln.h'^ = 0.146 ± 0.01. 

To fit SN at z < 2, DGP needs low matter density 
around ~ 0.2 which is 60% lower than best fit flm 
of DE. But with such a low we are not able to get 
the longer D{z\ss) measured by WMAP. It leads Ax^ ^ 5 
worse fit for DGP models in comparison with DE. Flat 
DGP models fail to the combined test of SN and CMB 
at around the 2 — cr confidence level. 

Though flat DGP type MG model fails to pass the 
test of SN and CMB, relaxing the fiat constraint helps 
DGP to fit dataset [13. By adding curvature, we are 
able to elongate D{ziss) with keeping low f2„j. If we put 
a few % curvature density to the total energy density, we 
have well-matched D{ziss)- For example, the curvature 
with fife = 0.05 gives DGP almost identical ~ 178 
of DE. But we do need a higher Hq for DGP around 
80km/s/Mpc which is 1 — ct level off from the Hubble con- 
stant measurement reported by HST 16] . It is also pos- 
sible that any other independent measurement of fi^ 
or Hq can also discriminate DE and DGP in the future. 

Inspired by DGP brane world model, the DGP type 
MG models can be extended generally I18!| . The modified 
Friedman equation can be written in the following ways 
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where the crossover distance Tc is defined as the ratio of 
5-dimensional to 4-dimensional Planck mass scales 
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The crossover distance Tc is a free parameter of DGP. If 
To is close to the current horizon scale, the acceleration 
of cosmic expansion is replicated without dark energy. If 
n = 2, we have known DGP theory and if n = cxd, then we 
have the cosmological constant. We find that DGP type 
MG models with n>2 shows good fit to SN+CMB. GEF 
test alone shows that the current GEF measurement does 
not rule out any theory. We start our estimation based 
upon the expected failure of test by GEF alone. 
Deviation parameters The combined test of GEF and 
GRF leads to the best opportunity to test standard grav- 
ity at large scales. The GEF and GRF can be simul- 
taneously determined at the weak gravitational lensing 
maps distorting the shape of galaxies systematically. The 
transverse shear 74 of the shapes of galaxies is gravita- 
tionally lensed by the intervening mass. The statistically 
measurable quantity 74 is generated by the GEF projec- 
tion of the transverse displacement driven by GRF. 

We group galaxies in n;, redshift bins and define the 
tangential shear 7^ in given bin i. The two point func- 
tion of this shear, {'ji^j), includes the information of both 
GEF and GRF. The power spectra in the spherical har- 
monic space Ci are used in our statistical estimation, de- 
fined by (7f™7j'"'') = Ci6wS„,,^,. We have Wnh + l)/2 
shear power spectra which are given by 0, |23, 0| 

= ^'^y drrWHr)W\r)A'i{k,z{r)) (3) 

where W(fi^ r) is the weight function |21| and the dinien- 
sionless power spectrum A_ is two point function of the 
quantity related to the net tangential deflection. 

The net tangential displacement by the gravitational 
lensing is given by the spatial gradient of the difference 
curvature perturbations = (<& — ^)/2, which is writ- 
ten as in real space (without (r) , $ and \1/ are in Fourier 
space.) 

Vj^$_(r) = V±[$(r) - *(r)]/2 (4) 

In the standard gravity frame, the difference curvature 
perturbation $_ is equal to $ with the condition ~ 
($ -I- ^')/2 = of ignorable anisotropy stress. The power 
spectrum A_ is given by 

A- = ;^(<f-<f->. (5) 

is not different from A$ in standard gravity. But in gen- 
eral MG models, A_ is not always equal to A<i,. What we 
measure by the weak lensing is the difference curvature 
perturbations 

We need three steps to calculate the measurable quan- 
tity $_: 1) to derive the matter perturbations, 2) to 
relate the matter perturbations to the curvature pertur- 
bations through non-dynamical constraint equation given 
by 00-component of Einstein equations and 3) to find the 
anisotropy sources which break the condition of (f>_(- = 0. 



We derive the matter perturbations by solving the fol- 
lowing differential equations, 

L + T^^m = -fc'*, (6) 

where the evolution of is sourced by the Newtonian 
constant \E'. Next the curvature perturbations <& is re- 
lated to ^iji by Poisson equation given by 

k^<l> = 4:^TGa^p„^S„^. (7) 

Then we get from the given condition $+. 

In DE model (standard gravity case) , there is no mod- 
ification to the basic form of Eq. and Eq. The ex- 
pansion history H(z) measured by GEF completely fixes 
In MG model, the modified 'J affects on the evo- 
lution of (5ni differently from standard gravity. Also the 
relation between <I> and is altered by the modified 
Poisson equation. We would get the distinct <!>_ of MG 
through possible modifications on Eq. (BJ Eq. [T] and the 
modified condition $4.. 

The general behavior of the deviation of $_ evolution 
from standard gravity can be deduced from the DGP 
type modification. In the DGP type model, the metric 
perturbations are modified in scalar-tensor type 

kH =- Wg[l + /3D^^(i/)]p„,^„, (8) 

Though we do not have the full solution of the 5 dimen- 
sional perturbation theory, we solve the 4D gauge the- 
ory on the brane based on the ansatz which closes the 
4D Bianchi identity. The exact form of /3^'^^{H) and 
P^'^^{H) depends on our ansatz on the 4D brane. 

In 4D DGP, we assume no significant contribution from 
the 5th dimensional projected Weyl fluid ini. This 
ansatz leads to no anisotropy stress on the brane, i.e. 
the condition = still holds 

f3tP{H)^P%^{H). (9) 

The reduced fundamental mass scale results in the in- 
crement of Gn variation which makes both $ and ^ 
increase. Shown in Fig ^ (upper long dash line in the 
bottom panel), $_ of 4D gets peculiar evolution at later 
time until $_ starts to decrease by the decay of Sm- This 
ansatz does not achieve the self-consistency. To neglect 
5th dimensional effects causes the unstable structure for- 
mation on the brane j^. But we still use this case to 
show the generality of our parameterization as a possible 
interesting case of MG which can boost higher. 

In static DGP, the time derivative to the spatial deriva- 
tive is ignored, which is a sound assumption for cos- 
mological structure formation. Then the non-negligible 
anisotropy stress tensor should be introduced 12]. The 
modification on the geometrical perturbations looks dif- 
ferent from the 4D DGP case as 

ps^aUcf^jj-^ = -/3^*°*''=(i/). (10) 
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While 4>+ = condition is broken by the significant con- 
tribution of the projected Weyl tensor, the deflected hght 
path by the gravitational lensing is unmodified. Shown in 
FigH (lower long dash line in the bottom panel), <!>_ gets 
the extra decrement in comparison with standard gravity 
case. Though <i>_ experience no modification, the weaker 

sourcing matter perturbation leads to decay of <i>_. 

Though DGP is found to have trouble in providing 
the stable structure formation, both opposite MG models 
help us to guess the general behavior of scalar-tensor type 
modification on As far as the scale dependence of 
GRF can be separated from the time dependence of GRF, 
we can design new deviation parameters by considering 
the modification effects on $+ and <i>_ like the DGP type. 

It leads us to consider the following empirical param- 
eterization as 



$' = (l + ZJ^a"")* 
*' = {1 + p^.a''")-^ 



(11) 



where /3$ and denote the constant strength of the 
deviation and Up determines the moment of the deviation 
from standard gravity. The shifted metric perturbations 

and well reconstruct the actual shape of observable 
$_ of each theory. 

The power rip needs to be positive to satisfy with the 
request of negligible modification effect at the matter 
dominated epoch. There is significant degeneracy to gen- 
erate among deviation parameters. As far as n^ is 
set to be reasonable value between 0.5 < Hp < 2.5, we 
are able to find a different pair of /3$ and to provide 
the approximately identical in the end. We choose 
np = 1.5 as a fixed constant and float /3$ and param- 
eters. 

We produce <i>_ based upon the real MG theory de- 
scribed above. Though we use the nearly identical input 
GEF of each theory, the difference curvature perturba- 
tions (f>_ of each theory are noticeably different from each 
other (the bottom panel of Fig . There are clear dis- 
tinct between 4D DGP and static DGP (upper long dash 
and lower long dash lines respectively) . The index n can 
denote the type of newly introduced non-linear energy 
momentum tensor; quadrature for n = 2 (DGP), cubic 
for n = 3 and quartic for n = 4. The bottom panel shows 
that as we increase n we would have closer of LCDM. 
We get $_ of n = 3,4 by using static ansatz only 
(dash and dash-dotted fines respectively). 

We try to reconstruct modified <&_ with (3 parameters 
without assuming any peculiar MG theory to begin with. 
The thin solid lines in the bottom panel of Fig ^ present 
$_ provided by /3 parameters. We find that our deviation 
parameters show great fit to any known MG theory. 
Forecast We float all known cosmological parameters 
without any prior information. Those are simultaneously 
determined by shear correlation functions and CMB 
power spectra (Planck). We take the same set of QCDM 
cosmological parameters in plj and add our /3 parame- 
ters with fixed np = 1.5. We consider a reference shear 
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FIG. 1: Top panel: the expected deviation parameter errors of /3$ 
and /3g, calculated by G2it weak lensing survey at fixed np = 1.5. 
Bottom panel: we show the observable structure formation by 
weak lensing <&_. Dotted line-standard gravity; two thick solid 
lines surrounding the standard gravity solution-translating 1 — a 
off curves of the top panel in space; two long dash lines-n = 2 
case, upper one for 4D DGP and lower one for static DGP; a dash 
lines and a dot dash line-n = 3 and n = 4 cases. We show the 
reconstructed 3>_ by our beta parameters at fixed np = 1.5. Four 
thin solid lines denote the reconstructed <I>— which show great fit 
to actual theories. 



survey, ground-based (G27r [23). The ground-based sur- 
veys can cover more sky. G27r represents the ground- 
based survey with a half-sky coverage. 

We calculate the uncertainties of parameters by using 
the first order Taylor expansion of the parameter depen- 
dence of two-point correlation functions. If the assump- 
tion of this linear response is available, the inverse of 
the Fisher matrix leads to the expected parameter errors 
with the given expected experimental errors on the power 
spectra. The contribution to the Fisher matrix from the 
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shear-shear correlations is given by 

^= E ^c;,^^^wf-^^c;;^;^<^^\ (12) 

i,il,i2,i3,i4 

where the subscript ,p denotes differentiation with re- 
spect to parameter Op. We use W to denote the inverse 
of the total covariance matrix: 

>V=(S + N)-i. (13) 

Note that W is an easily invertible block diagonal matrix: 
Wirm,i'm'i' = Smm' Sw . Thc total Fishcr matrix is 
given by summing F'^'^ with the Fisher matrix for the 
unlensed CMB data, F^^b 

We take the standard gravity model (DE) as our fidu- 
cial one. If DE is a righteous modeling for cosmic accel- 
eration, then we would be interested in the upper bound 
of standard gravity. On the other hand, if MG is a cor- 
rect one to explain cosmic acceleration, then we need to 
know if proper MG model can be discriminated from the 
standard case and from the other MG models. 

Standard gravity points (0,0) on the (3 space. We es- 
timate the expected errors with G2tt in the top panel of 
FigD The errors on f3 parameters indicate the fractional 
difference of modified gravitational constant from stan- 
dard Gat at the present time. We get the upper bound 
for effective Gn around 8.8% level of accuracy with G2tt 
at the present time (effective fractional error ac^ for ^- 
is defined by (erg • ■ (fpY^"^ /2.). The percentage er- 
ror for Gn ranges from 3% to 15% with 0.5 < n^ < 2.5. 
The estimation of Gn error gives us the rough idea of 
constraints on the standard gravity. But considering the 
phenomenological aspects of weak lensing, it would be 
better to translate the beta parameter constraints on the 
space where we can see the evolution of the structure 
formation. 

In the bottom panel of Fig ^ we translate the con- 
straints of beta parameters on $_ plane. The result is 
really encouraging. Approximately 10% uncertainty in 
Gn results in about 2% level of constraints on evolu- 
tion at later time. If the standard gravity is true model of 



our universe even at such a large scale, then we can pre- 
dict the structure formation within 2% level of accuracy 
by determining GEE precisely. It could get worse by all 
the systematic error expected with G27r, but it is recom- 
mended to combine all other GEE measurements other 
than weak lensing (SN and baryon oscillation), which will 
lead to much stronger constraints on <f>_ evolution after 
all. 



The meaning of constraints on $_ we get can be seen 
by comparing the estimated error with the level of depar- 
ture from the standard gravity predicted by DGP type 
MG model (the bottom panel of Fig ^ . Obviously both 
n = 2 cases are distinguishable from the standard gravity. 
Also even higher n cases are marginally discriminable. 
Considering all known MG models, we can conclude that 
the constraints on the standard gravity estimated by us- 
ing G27r are meaningful ones to test the standard gravity 
at large scales in the future. 

Conclusion We design the deviation parameters and cal- 
culate the expected errors of those parameters by using 
the combination of weak lensing and CMB. From the ex- 
pected errors on (3 parameters, we calculate the upper 
bound of effective Gn{^) around 10% accuracy at the 
present time by measuring deviation from standard 
gravity within 2% accuracy. Shown in Fig^ it is well be- 
low the level of a departure of $_ reported by any known 
MG model so far. DE would be strongly supported by 
this test with nearly zero [3 parameters measured. On 
the other hand, any measurement of non-zero (3 param- 
eters could lead us to shift to MG modeling from DE to 
explain cosmic acceleration, or to look for any other mi- 
nor broken standard gravity with keeping DE modeling 
as it is. 



Acknowledgments: We thank Christopher Gordon, 
Wayne Hu, Ignacy Sawicki and Xiaoming Wang for help- 
ful discussion. YSS is supported by the U.S. Dept. of 
Energy contract DE-FG02-90ER-40560. This work was 
carried out at the KICP under NSF PHY-0 114422. 



[1] S. Perlmutter et al. (Supernova Cosmology Project), As- 

trophys. J. 517, 565 (1999). 
[2] G. R. Dvali, G. Gabadadze, and M. Porrati, Phys. Lett. 

B485, 208 (2000). 
[3] C. Deffayet, Phys. Lett. B502, 199 (2001). 
[4] S. M. Carroll, V. Duvvuri, M. Trodden, and M. S. Turner, 

Phys. Rev. D70, 043528 (2004). 
[5] A. Lue, R. Scoccimarro, and G. D. Starkman, Phys. Rev. 

D69, 124015 (2004). 
[6] Y.-S. Song, Phys. Rev. D71, 024026 (2005). 
[7] C. Schimd, J.-P. Uzan, and A. Riazuelo, Phys. Rev. D71, 

083512 (2005). 

[8] L. Knox, Y.-S. Song, and J. A. Tyson (2005), astro- 



ph/0503644. 

[9] M. Ishak, A. Upadhye, and D. N. Spergel (2005), astro- 
ph/0507184. 

[10] E. V. Linder, Phys. Rev. D72, 043529 (2005). 
[11] I. Sawicki and S. M. Carroll (2005), astro-ph/0510364. 
[12] K. Koyama and R. Maartens, JCAP 0601, 016 (2006). 
[13] A. G. Riess et al. (Supernova Search Team), Astrophys. 

J. 607, 665 (2004). 
[14] D. N. Spergel et al. (WMAP), Astrophys. J. Suppl. 148, 

175 (2003). 

[15] M. Fairbairn and A. Goobar (2005), astro-ph/0511029. 
[16] W. L. Freedman et al., Astrophys. J. 553, 47 (2001). 
[17] G. Bernstein (2005), astro-ph/0503276. 



5 



[18] G. Dvali and M. S. Turner (2003), astro-ph/0301510. 
[19] W. Hu, ApJ 522, L21 (1999). 

[20] M. Takada and B. Jain, Mon. Not. Roy. Astron. Soc. 

348, 897 (2004). 
[21] Y.-S. Song and L. Knox, Phys. Rev. D68, 043518 (2003). 



[22] K. Koyama (2006), astro-ph/0601220. 
[23] J. A. Tyson, D. M. Wittrnan, J. F. Hcnnawi, and D. N. 
Spergel, Nucl. Phys. Proc. Suppl. 124, 21 (2003). 



